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Abstract. On 18 November 1993, we observed the QSO 
1150+812 and the BL Lac object 1803+784, nearly 
15° apart on the sky, with a VLBI array, simultaneously 
recording at 8.4 and 2.3 GHz. Using difference astrometry 
techniques, we determined the coordinates of 1803+784 
relative to those of 1150+812 to be 

Aa = 6 fe 7 m 33U8469 ± K00020 
AS = -2°30'25'.'13557±0'.'00075 



These standard errors contain estimated contributions 
from the propagation medium, the effects of source opac- 
ity, and the possible misidentification of the reference 
point in some of the maps. We obtained comparably accu- 
rate and consistent relative positions of the radio sources 
with GPS-based ionosphere corrections, thus demonstrat- 
ing that dual-frequency observations are not required 
for state-of-the-art accuracy in VLBI astrometry. These 
results also demonstrate the feasibility of using single- 
frequency measurements for sources separated by 15° on 
the sky, and open the avenue for the application of this 
technique on a full-sky scale. 



Key words: astrometry - techniques: interferometric 
- quasars: individual: 1150+812 - BL Lac: individual: 
1803+784 



1. Introduction 

The establishment of a radio reference frame with sub- 
milliarcsecond accuracy has been a major goal of as- 
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trometrists for the last several decades. Centimeter- 
wavelength very-long-baseline interferometry (VLBI) 
group-delay astrometry of extragalactic radio sources rou- 
tinely provides precisions at the milliarcsecond (mas) 
level, thus allowing a celestial reference frame to be built 
with corresponding accuracy (e.g., Ma et al. 1998[ ). The 
use of phase-delay difference astrometry (see, e.g., Shapiro 
et al. 1979] ) allows us to determine angular separations 
between pairs of radio sources at the submilliarcsecond 
level. One important advantage of phase-delay astrome- 
try is that we can identify with sufficient accuracy suitable 
reference points within the structures of the radio sources 
whose relative positions we wish to determine, whereas in 
group-delay astrometry suitable reference points are not 
so easily identified from epoch to epoch with the desired 
accuracy. 

With the phase-delay technique, the differences be- 
tween phase delays for two radio sources are used to de- 
termine their relative positions, because these differences 
have reduced sensitivity to unmodeled effects. For sources 
nearby to one another on the sky, this technique yields 
statistical standard errors of only a few microarcseconds 
(/ias), as in the case of quasars 1038 +528 A and B (Mar- 
caide & Shapiro 1983 ; Marcaide et al. 1994 ), whose compo- 
nents are separated by only 33'.' However, the overall stan- 
dard errors are dominated by inaccurac ies in the reference- 
point identification (e.g., Rioja et al. 1997 ). For sources 
with larger separations, the main contributions to the as- 
trometric standard error in the relative position of two 
sources comes from uncertainties in the coordinates of the 
source chosen as the reference, in the value used for UT1- 
UTC, and in the effects of the propagation medium. For 
radio sources with separations ranging from 0.5° to 7° on 
the sky, standard errors in relative positions of about 0.1- 



0.3 mas have been obtained (Bartel et al. 1986; Guirado 
et al. |l995ajb, |1998|; Lara et al. |1996| Ros et al. [1999^. 



A limiting factor in centimeter-wavelength VLBI as- 
trometry is the uncertainty in the contribution of the iono- 
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sphere to the astrometric observables, even though its dis- 
persive character makes this contribution scale as v~ 2 and, 
in principle, allows it to be determined accurately (see, 
e.g., Thompson et al. 1986 ). One strategy to take advan- 
tage of this scaling is to make simultaneous VLBI measure- 
ments in two frequency bands. The main disadvantage of 
this option is that the (fixed) bandwidth of the recording 
equipment has to be split between two frequency bands, 
decreasing the signal-to-noise-ratio (SNR) for each. Alter- 
natively, one may compute corrections based on estimates 
of the ionosphere total electron content (TEC) obtained 
independently from Faraday-rotation measurements and, 
more recently, from the Global Positioning System (GPS). 
The advantage of the latter approach (e.g., Guirado et al. 
1995b| , Ros et al. |2000| ), is that only single-band VLBI 
observations are needed, avoiding the loss of sensitivity 
mentioned above. 

In this paper, we study the applicability of the phase- 
difference technique to the strong radio sources 1150+812 
and 1803+784, separated on the sky by almost 15°. We 
show that reliable phase connection is feasible at such a 
large angular separation, and estimate the relative posi- 
tion of the two sources with submilliarcsecond accuracy. 
We compare the estimates of the relative angular separa- 
tion that result from use of two different methods of re- 
moving the ionosphere contribution based on two different 
types of data, namely, phase delays from dual-frequency- 
band VLBI measurements and TEC values from GPS 
measurements. The estimates from the two methods agree 
to within the standard errors from each method, showing 
that single-frequency astrometric VLBI experiments can 
be confidently carried out. 

In Sect. |^, we briefly describe the observations; in 
Sect. || we describe the radio sources, and in Sect. f| the 
data reduction process. In Sect. || we discuss our estimate 
of the relative position of the source 1803+784 with re- 
spect to 1150+812, and in Sect. ^ we carry out a sensitivity 
analysis of the estimated position of 1803+784 to errors in 
other model parameters. Finally, in Sect. [7] we summarize 
our main results and discuss their implications. 
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Fig. 1. Hybrid maps of 1803+784 at 8.4 (top) 
and 2.3 GHz (middle and bottom). Contours are 
0.5,1,2,4,8,16,32,64, and 90% of the peaks of brightness, 
which are 1.48 Jy/beam and 1.35 Jy/beam for 8.4 and 
2.3 GHz, respectively. The sizes of the restoring beams, 
shown shaded at the bottom left corners, are 0.52 x 0.49 
mas (PA=61°) at 8.4 GHz and 2.21 x 1.68 mas (PA=29°) 
at 2.3 GHz. East is left, and north is up. The identifica- 
tion of the components in the 2.3 GHz map is based on 
an overresolved map (not shown here). 



2. Observations 

We observed the radio sources 1150+812 and 1803+784 
on 18 November 1993 for 12 hours, in right circular polar- 
ization, simultaneously recording at two frequency bands 
(X-band w 8.4GHz and S-band w 2.3 GHz). We used the 
following radio telescopes (in parenthesis we give the site 
symbol used in this paper, the diameter, and the location 
of the telecope): Effelsberg (B, 100m, Germany); Medic- 
ina (L, 32m, Italy); Onsala (S, 20m, Sweden); Fort Davis 
(F, 25m, Texas); Hancock (H, 25m, New Hampshire); 
North Liberty (I, 25m, Iowa); Owens Valley (O, 25m, 
California); Los Alamos (X, 25m, New Mexico); and the 
phased VLA (Y, 130 m equivalent, New Mexico). The Eu- 
ropean antennas (B, L, and S) recorded in Mark III mode 



A, covering a total bandwidth of 56 MHz, via seven ad- 
jacent channels spanning 8,403 to 8,431 MHz and seven 
such channels spanning 2,273 to 2,301 MHz. The Ameri- 
can antennas used the Mark III recording system in mode 

B, covering a total bandwidth of 28 MHz, with four con- 
tiguous channels spanning 8,403 to 8,419 MHz and three 
such channels spanning 2,289 to 2,301 MHz. We used an 
observing cycle consisting of 2 min observing 1150+812, 
1.5 min antenna slewing, 2 min observing 1803+784, and 
1.5 min antenna slewing back to 1150+812, making a total 
cycle duration of 7 min. 

The data were correlated at the MPIfR correlator in 
Bonn, Germany. Since both sources were strong (> 1 Jy) 
at the epoch of observation, we detected each of them 
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Fig. 2. Hybrid maps of 1150+812 at 8.4 (left) and 2.3 GHz (right). Contours are 0.5,1,2,4,8,16,32,64 and 90% of the 
peaks of brightness, which are 0.90 Jy/beam and 0.74 Jy/beam for 8.4 and 2.3 GHz, respectively. The sizes of the 
restoring beams, shown shaded at the upper right corners, are 0.61 x 0.55 mas (PA=15°) at 8.4 GHz and 2.33 x 1.67 
mas (PA=-56°) at 2.3 GHz. The identification of the components in the 2.3 GHz map is based on an overresolved map 
(not shown here; see text). 



with high SNR within the 2 min integration time, at both 
frequencies and for each baseline, except those baselines 
involving the Onsala antenna. After correlation, we ex- 
ported the VLBI observables (fringe amplitude and phase, 
group delay, and phase-delay rate) for the reference fre- 
quencies fx = 8413 MHz and us = 2295 MHz, and cal- 
ibrated the fringe amplitudes using the information pro- 
vided by the staffs at the observing antennas. In our astro- 
metric analysis, we discarded the data from Onsala due to 
the low SNR of the detections with this antenna; and from 
the phased VLA, since this antenna recorded only at 8.4 
GHz and was used solely to increase the SNR of our hy- 
brid maps. In addition, we discarded data from Effelsberg 
for the interval 18:30 until 23:28, since this antenna, too, 
only observed at 8.4 GHz during that interval. To obtain 
our hybrid maps, we self-calibrated, Fourier inverted, and 
CLEANed the visibility data using the Caltech program 



DIFMAP (Shepherd et al. 1995) 



3. The Radio Sources 

The radio source 1803+784 has been optically ide ntifie d 
with a BL Lac object at z=0.684 (Stickel et al. 



1991|) 



This source has a nearly flat spectrum from millimeter to 
decimeter wavelengths (Strom & Biermann 1991) and dis- 
plays many interesting properties (Biermann et al. 1981), 



including rapid variability, compactness, and strong X-ray 
emission. 



Our 8.4 and 2.3 GHz maps of 1803+784 (Fig. 0) show a 
jet pointing westward. The hybrid map at 8.4 GHz shows 
the jet extending about 3 mas, with three components, 
labeled XI, X2, and X3. Most of the emission comes from 
components XI, identified as the core, and X2, about 1.5 
mas west of XI. Component X3 lies about 2.8 mas west 
of XI and is much fainter than XI and X2. Our map at 
2.3 G Hz als o shows the jet, w hich, as in the maps of Ros 
et al. ( 1999] ) and Eckart et al. (1987), extends up to about 
25 mas and is slightly bent towards the southwest. At 2.3 
GHz, most of the emission is concentrated within 3 mas of 
the peak, but emission is also discernible at points about 
5, 8, and 10 mas west of SI (Fig. 0, middle). These lat- 
ter bright regions, especially the innermost and the outer- 
most, may correspond to real features of the jet, but like 
the questionable features farther out, they are not rele- 
vant for our astrometric purposes and, therefore, are not 
labeled. 

The radio source 1150+812 is a QSO with redshift z = 
1.25 (Hewitt & Burbidg e jL993D , For this source, our 8.4 
and 2.3 GHz maps (Fig. |2|) indicate a curved jet directed 
to the south and west. The 8.4 GHz map contains at least 
four components, labeled XI through X4, from north to 
south. Most of the emission comes from the components 
XI and X2, which are separated by ~ 1.5 mas. Extended 
emission is still discernible at a distance of ~ 4 mas from 
XI (components X3 and X4). At 2.3 GHz, the emission is 
concentrated within the innermost 5 mas, but emission is 
also evident farther to the southeast. 
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Our maps of 1803+784 and 1150+812 are in good 
agreement with maps of similar resolution from 1995 pub- 
lished by Fey & Chariot ( |1997| ). 



4. Astrometric Data Reduction 

Our strategy in the astrometric data analysis was to ob- 
tain a set of "connected" phase delays that could be an- 
alyzed via weighted-least-squares to estimate the relative 
position of the two sources. 

The observables in our experiment, after correlating 
and fringe-fitting the data, were the interferometric phase, 
4>, the group delay, tq, and the phase-delay rate, f^. We 
based our analysis primarily on 0, converted first to the 
phase delay, r^, via removal of the "27rn" ambiguities 
(where n is an integer) . The delay equivalent of 2 n in 
phase is ~ 120 ps at 8.4 GHz and ~ 440 ps at 2.3 GHz. 
To remove these ambiguities we used a "phase connec- 
tion" technique (see, e.g., Shapiro et al. 1979) that can 
be outlined as follows: We constructed a model of the 
phase delays, based on a model of the geometry of our 
interferometric array and the radio sources (Table |l|); the 
propagation medium; and the clock behavior at each sta- 
tion relative to a reference. We estimated the parameters 
of this model via weighted-least-squares analysis of the 
observed phase-delay rates and group delays. We us ed an 
improved version of the program VLBI3 (Robertson 1975) 
to carry out this weighted-least-squares analysis. The re- 
sultant model of the phase delays was accurate enough to 
allow phase connection, i.e., the elimination of the 2-nn 
ambiguities for almost all of the observations (see below), 
via a suitable iterative scheme that took advantage of 
the closure relations over triangles of baselines. Since we 
phase-connected the data independently for each source, 
we used weighted-least-squares to verify - successfully - 
the consistency of the "overall resolution of a mbigui ties" 
of the delays for the two sources (see Ros et al. 1999 for a 
detailed discussion of this step). During our analysis, we 
discarded all data from the Hancock antenna, since our 
phase connection failed at both frequencies. We also dis- 
carded the data from the Medicina antenna from 18:30 
to 23:48, since in this time period our phase-connection 
appeared to be unreliable. (These failures are likely the 
result of some combination of tropospheric, ionospheric, 
and instrumental effects.) 

1.1. Source- Structure Contribution 



The extended structure of even a "compact" radio source 
often makes a significant contribution to the phase de- 
lays. Such a contribution depends on the point chosen as 
a reference on the map. Identifying a reference point in a 
reliably epoch-independent manner is crucial for the use 
of our method to compare positions obtained from dif- 
ferent epochs. For each source and frequency band, our 
procedure was to choose the peak of brightness (POB) as 



the reference point. To find the POB, we constructed fine- 
grained maps of the two radio sources using pixel size 0.01 
mas at 8.4 GHz and 0.03 mas at 2.3 GHz. We identified the 
brightest pixel in each such map as the POB and defined 
a new coordinate system with that reference point as its 
origin. We then computed from these maps the structure- 
phase contribution to each phase delay. These contribu- 
tions - up to ~ 25 ps at 8.4 GHz and ~ 55 ps at 2.3 GHz 
- were removed from the phase delays, to effect a point-like 
source at each POB. The standard deviations assigned to 
each POB location included uncertainties due to (a) the 
pixel size used for each map, which was about 0.01 (0.03) 
mas at 8.4 (2.3) GHz, for each source; and (b) SNR for the 
peak, which was about 0.04 (0.27) mas at 8.4 (2.3) GHz, 
for each source. Note, however, that when the delays at 
the two frequencies are combined to form "plasma-free" 
delays, errors in the phase delays at 2.3 GHz are scaled 
down by a factor [{v x /v s ) 2 - 1]~ ~ 0.08 (see Sect. O), 



and therefore the uncertainty of the POB at 2.3 GHz is 
less significant than that of the POB at 8.4 GHz. The root- 
sum-squares of these contributions, at each frequency, are 
indicated as the first (8.4 GHz) and second (2.3 GHz) en- 
tries in Table and are dominated by the relatively low 
SNR of the hybrid maps. 



4-2. Neutral Atmosphere Contribution 

The neutral atmosphere (primarily the troposphere) adds 
an extra delay to the incoming radio waves, the equiva- 
lent of up to a few meters in pathlength. We monitored 
the pressure, relative humidity, and temperature at each 
observing antenna to track the atmosphere behavior dur- 
ing the observations. We used the model by Saastamoinen 



1993, in which the atmosphere is separated into two com- 



ponents: a dry component and a wet component (due to 
the water vapor in the atmosphere). Based on this model, 
we calculated a priori values for the delays for the wet 
and dry atmosphere components in the zenith direction 
f or ea ch antenna site (Table [l]) , and then used the Chao 
( 1974 ) mapping function to determine delays at other el- 
evations. This mapping function agrees to within about 
1 cm with ray-tracing computations (Davis et al. 1985 ) 
for antenna elevations larger than 20°; all of our obser- 
vations had elevation angles above this limit. To specify 
adjustments to the combined (dry and wet) atmosphere 
delays during the observations, we used a piecewise- linear 
function characterized by zenith delay values at epochs 
about three hours apart. Errors in the combined neutral- 
atmosphere delays - mainly due to the wet troposphere - 
are likely to be up to 0.1 ns/hr (D. Lebach, priv. comm.). 
Since the wet atmosphere zenith delay fluctuates approx- 
imately as a random walk (Treuhaft & Lanyi 1987| ), an 
error of 0.1 ns/hr transforms into a standard error of 
V3 x 0.1 w 0.17 ns every three hours. Therefore, in our 
sensitivity analysis, we allowed the atmosphere parame- 
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Table 1. A priori values of the parameters of the theoretical model used in the estimates of the source positions 



Source coordinates (J2000.0) a 



1150H 
1803H 



812 
-784 b 



a 
a 



ll n 53 m 125499188 
18 h 00 m 45f683922 



: 80°58'29"15453 
: 78°28'04"01852 



Antenna site coordinates and atmosphere parameters 



Radio 

telescopes 3 



B 
F 
I 
L 

O 

X 



Effelsberg 
Fort Davis 
North Liberty 
Medicina 
Owens Valley 
Los Alamos 



Cartesian coordinates [ml 



X 

4033947.51 
-1324009.08 
-130872.19 
4461370.05 
-2409150.05 
-1449752.31 



Y 

486990.46 
-5332181.98 
-4762317.13 

919596.77 
-4478573.27 
-4975298.60 



4900430.75 
3231962.49 
4226851.05 
4449559.16 
3838617.42 
3709123.95 



Axis 



Mean atmosphere 



offset [m] zenith delay c [ns] 



0.00 
2.12 
2.13 
1.83 
2.13 
2.14 



7"dry 

7.5 
6.4 
7.5 
7.7 
6.8 
6.2 



7~wct 

0.1 
0.1 
0.1 
0.2 
0.1 
0.1 



Earth Tides 

Radial Love number d , h = 

Horizontal Love number d , 



: 0.60967 
1 = 0.085 



Precession Constant (J2000.0) c 

p=5029".0966/Julian Century 
Mean obliquity (J2000.0) c 

e =23°26'2l"448=8438l"448 



Tidal lag angle, 6 = 0°.0 



Polar Motion and UTl 

Interpolate d from daily values 
from IERS (199C) 



Nutation 



Interpolatio n of] ATT 1980 Nutation Model 
3q) daily corrections 



with IERS ( 



1996 



We took 0.058 ms and 0. 140 m as to be the standard errors in the right ascension and declination, respectively, of 1150+812, 
twice the quoted IERS (1996) errors. We also used as standard errors 2 cm in each coordinate of every antenna location 
(conservatively chosen to be at least twofold larger than the errors given by IERS), 0.7 mas in each pole coordinate, and 
0.04 ms in UT1-UTC. We calculated the location of each antenna by linear extrapolation, using the velocity given by the 
IERS, and applied to the 1993.0 position, also given by the IERS. 

The coordinates of 1803+784 serve only as an initial estimate, and do not affect our results. 

Mean values of sampled zenith delays, computed using the meteorological values provided by the staff at the observing 
antennas, 



Dahlen ( J1976 V 

Lieske et aT~(1977j) 



ters at each antenna location to vary with this standard 



4-3. Ionosphere Contribution 

As explained in Sect. [j], we removed most of the iono- 
sphere contribution in two alternative ways, first by us- 
ing our dual-frequency-band observations and, second, by 
using the total electron content (TEC) values d educe d 
from GPS measurements (see, e.g., Sardon et al. 1994 ). 
For the latter removal, we followed the same procedure as 
described in Ros et al. (2000), modeling the ionosphere 
as a thin layer located at an altitude of 350 km. We used 
GPS data from Wettzell (Germany) to obtain TEC values 
for Effelsberg and Medicina, and from Goldstone (Califor- 
nia) to obtain TEC values for Owens Valley, Los Alamos, 
and Fort Davis; we were fortunate to have a GPS antenna 
collocated with our VLBI antenna at North Liberty to 
obtain TEC values for this site. In Fig. || we compare 
the estimated ionosphere delays as deduced from both 
dual-frequency-band observations (8.4 GHz/2.3 GHz de- 
lays) and GPS-based measurements (TEC delays), for a 
representative subset of baselines, intracontinental (upper 
two plots) and intercontinental (lower two plots). Dual- 
frequency corrections refer only to baselines, not individ- 
ual antennas. On the other hand, TEC ionosphere correc- 



tions are calculated for individual antennas. Therefore, for 
each antenna there can be a constant offset between the 
TEC ionosphere correction and the dual-frequency one as 
a result of instrumental effects, but, for clarity, we sub- 
tracted the mean difference between each of the two cor- 
rections (Effelsberg, 1.35 ns; Medicina, 0.55 ns; North Lib- 
erty, 7.80 ns; Fort Davis and Los Alamos had no offset). 
We assumed a statistical standard error for the vertical 
TEC at each GPS antenna of 2 x 10 16 mT 2 (Sardon, priv. 
comm.), corresponding to a standard error of ~40 ps for 
the (vertical) phase delays at 8.4 GHz. For each antenna, 
this error is multiplied by the value from the mapping 
function for the elevation angle of each radio source (here, 
the secant of the zenith angle at the ionospheric point) 
of each radio source. The resultant mapped TEC-based 
ionosphere corrections had standard errors ranging from 
~70 ps to ~ 120 ps, highly correlated from point to point 
(Fig. |). The standard errors for our dual-frequency-band 
corrections are the appropriate combination of the sta- 
tistical errors of the phase delays at each frequency, and 
ranged from ~ 5 ps for the Effelsberg-Medicina baseline 
up to ~ 30 ps for the Fort Davis-Owens Valley baseline. 
The root-mean-square (rms) of the differences between the 
corrections obtained by the two methods ranged from ~ 20 
ps up to ~ 90 ps, depending on baseline. The maximum 
difference was for intercontinental baselines and was ~ 110 
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ps, when it was night in Europe and noon in North Amer- 
ica. The level of agreement in the results from the two 
methods is gratifying, and foreshadows the agreement in 
the astrometric results obtained below with the two meth- 
ods. 



4-4- Opacity Contribution 

When we use phase delays at both 8.4 and 2.3 GHz (see 
Sect. 4.3), the reference points on the 2.3 GHz maps should 
be chosen to correspond to the same points in the sky as 
those on the 8.4 GHz maps. Unfortunately, opacity ef- 
fects may introduce an offset between the POB for one 
frequency band with respect to the POB for the other fre- 
quency band (see, e.g., Marcaide & Shapiro 1984). Since 



In this dual-frequency-band method to reduce the 
ionosphere contribution, the resultant "ionosphere-free" 
phase delays have the form: 



T <t>Jr 



R-l 



jet components are less self-absorbed than core compo- 
nents, their associated peaks should be less affected by 
opacity effects and should therefore be better suited for 
use as guides for the alignment of the maps at 8.4 and 
2.3 GHz. We registered the peaks of the less self-absorbed 
components in our 8.4 GHz maps with their corresponding 
peaks in (twofold overresolved) 2.3 GHz maps. 

For 1803+784, we identified component X3 at 8.4 GHz 
with component S3 at 2.3 GHz (Fig. [j]) on a twofold over- 
resolved map (not shown). Based on this map, we con- 
cluded that blended components SI, S2, and S3 at 2.3 
GHz correspond to components XI, X2, and X3, respec- 
tively, in the 8.4 GHz map. This identification implied 
that in our maps the POB of 1803+784 at 2.3 GHz (SI) 
is offset by ~ 0.6 mas westward from the POB at 8.4 GHz 
(XI). For 1150+812, we also used a twofold overresolved 
map at 2.3 GHz to make a plausible registration of the 
maps for the two bands. We identified components X3 
and X4 at 8.4 GHz with components S3 and S4 at 2.3 
GHz (Fig. |), implying a shift of the POB at 2.3 GHz 
(SI) of ~ 0.5 mas northward from the peak of the core 
in the 8.4 GHz map (XI). To estimate the standard error 
in our determination of those shifts we followed a similar 
procedure to that described for the determination of the 
standard error associated with the reference-point iden- 
tification (POB), but applied to components X3 and S3, 
for each source. This procedure yielded a standard error 
of 0.4 mas for each source. Finally, to take into account 
the possible additional shifts due to opacity effects of the 
POB of component S3 with respect to X3, for each source, 
we increased the estimated standard error in each coordi- 
nate of the registration estimate from 0.4 mas to 0.7 mas. 
Estimates of the shift between the POBs at 8.4 and 2.3 
GHz (where opacity effects are expected to be larger than 
in the jet components) of oth er radio sources - 103 8+52 8 
A (Marcaide & Shapi ro |l984| ), 1226+023 (Chariot [1993] ), 
1901+319 (Lara et al. 1994) - are also no greater than 0.7 
mas. Thus, our estimates of these shifts for 1803+784 and 
1150+812 seem reasonable. 



where the subscripts 1 and 2 refer to the two frequency 
bands, with V\ > V2 and R = {v\/v2) 2 - In our case, R ~ 
(8.4/2.3) 2 ~ 13; hence standard errors in the phase delays 
at 8.4 GHz are scaled up by a factor R/(R — 1) ~ 1.08, 
whereas corresponding errors in the phase delays at 2.3 
GHz are scaled down by l/(R— 1) ~ 0.08 in their effect on 
T0,frec- Likewise, uncertainties accounting for registration 
errors due to opacity effects have to be scaled by 1/(R— 1). 
Thus, the contribution to the total error budget due to the 
use of dual-frequency-band measurements, in particular 
the seemingly large error at 2.3 GHz discussed above, is 
not significant (see entry "Opacity effects" in Table ||) . 

Table 2. Contributions to the standard errors of the esti- 
mates of the relative coordinates of the POB at 8.4 GHz 
of 1150+812 and 1803+784 



Effect 


~ a 

(TAa 


OAS 




mas 


mas 


8.4 GHz source structure 13 


0.07 


0.07 


2.3 GHz source structure 


0.04 


0.04 


Opacity effects c 


0.08 


0.08 


Statistical errors 


0.41 


0.62 


Software limitations 


0.40 


0.40 


Root-sum-square 


0.58 


0.75 


Contributions to the standard errors of the estimates of the 



coordinates of 1803+784 minus those of 1150+812, namely 
relative right ascension, Aa, and relative declination, AS. 
3 Root-sum-square (rss) of the standard errors due to the use 
of data at 8.4 GHz. Each value has been scaled by 
R/{R - 1) « 1.08. See Sect. 

: rss of the standard errors due to the use of data at 2.3 GHz, 
which includes estimated errors in the "registration" of the 
2.3 GHz and 8.4 GHz maps, taken to be 0.7 mas for each 
sourc e. E ach value has been scaled by 1/(7? — 1) « 0.08. See 
Sect.fO. 

1 We estimated ctaq in milliseconds (ms). To transform this 
estimate into milliarcseconds (mas), we multiplied by 

15 • COS <5l8Q3+784 • 

; See Sect. |. 



5. Relative Position of the Two Sources 

We obtained a final set of phase delays for each source by 
correcting for source-structure, opacity-effect, and iono- 
sphere contributions as described above. We then formed 
a set of differenced phase delays by subtracting the delay 
for each observation of 1803+784 from the correspond- 
ing delay for the previous observation of 1150+812. The 
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Fig. 3. Ionosphere contribution, in nsec at 8.4 GHz, as obtained from our dual-frequency-band (8.4/2.3 GHz) mea- 
surements (filled circles) and from the GPS-based measurements (open circles), for a representative set of baselines 
including both intracontinental and intercontinental baselines. The standard errors shown are for the TEC-inferred 
delays, and are much larger than those from dual-frequcncy-band measurements. All four plots start at 00:00, since 
until approximately that time Effelsberg observed only at 8.4 GHz, and thus lacked dual-frequency-band corrections. 



use of differenced phase delays is in general more effective 
the closer together the sources are in the sky, since dif- 
ferencing results from neighboring observations tends to 
cancel effects which for each source alone cannot be accu- 
rately enough described by theoretical models. The best 
such pair of sources for such cancellation so far studied 
is 1038+528 A and B (Marcaide & Shapiro |1983|), be- 



cause this pair simultaneously lies inside the beam of each 
antenna, yielding almost complete cancellation of several 
sources of error. For sources separated by increasingly 
larger angular distances, the cancellation of unmodeled 
effects in general lessens, due directly to the increase in 
angular separation of the sources and indirectly to the in- 
crease of the cycle length of the observations. 

From a weighted-least-squares analysis of the dif- 
ferenced phase delays, we estimated the coordinates of 
1803+784 relative to those of 1150+812 (see Table |). In 
this analysis, we also included the (undifferenced) phase 
delays of 1150+812, suitably weighted, to estimate the 
relative behavior of the site clocks. (The opposite scheme, 
i.e., using the differenced phase delays of 1803+784 minus 
those for 1150+812, and the undifferenced phase delays of 
1803+784, did not alter our estimate significantly.) The 
total number of parameters whose values were estimated 
was 26, and include those of a third-degree polynomial to 
model the clock behavior at each station (except for North 



Liberty, which was taken as our reference station), the co- 
ordinates of 1803+784, and four atmosphere parameters 
(see Sect. 4.2) for the Effelsberg antenna, since the in- 
formation on the atmosphere conditions provided for this 
site was very sparse, and thus of doubtful utility. Nonethe- 
less, the estimated atmosphere parameters for Effelsberg 
agreed with the a priori values within 0.11 ns. 

Our choice of 1150+812 as the reference source was 
motivated by its being one of the defining sources of the 
quasi-inertial International Celestial Reference Frame. As 
a test of the effect of this choice, we solved for the coor- 
dinates of 1150+812-keeping those of 1803+784 fixed at 
their IERS values-and obtained different results for Aa 
and AS. As a simple check on the significance of those 
differences - one that avoids the effects of the reference 
system in which the coordinates are expressed - we cal- 
culated the arclength between the two sources for each of 
the two solutions. We expected to obtain nearly the same 
result for the two cases and, indeed, the difference was a 
mere 5 /ias, negligible compared to the overall standard 
errors of our results. 

For our nominal, dual-frequency-band solution de- 
scribed above, we scaled separately the standard devia- 
tions of the differenced phase delays and of the phase de- 
lays of 1803+784, so that for the data for each type, and 
for each baseline, the weighted-mean-square of the post- 
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Table 3. Phase-delay-based estimates of the coordinates of 1803+784 minus those of 1150+812 



Ionosphere correction method 



Aa~6 h 7 m 3S18473 c 
A5 + 2° 30'25'.'13601 c 



D ual- frequency 
(mas) 



GPS-derived TEC 
(mas) 



-0.12 + 0.58° 
-0.44 + 0.75 



+0.12 ±0.78 a 
-0.31+0.82 



a Estimates of the relative position of the two sources, minus the reference separation, from a weighted-least-squares analysis of 
the phase delays, after corrections for source-structure, source opacity, and ionosphere effects, the last via dual-frequency-band 
observations. The estimated standard errors include statistical standard errors and errors due to opacity effects, to registration 
of the maps, and to software limitations. See Table ^ and Sect. ^| 

b Estimates of the same relative position as above, except that only phase delays at 8.4 GHz were used and that the 
corrections for ionosphere effects were obtained from GPS-based estimates of TEC. 

c The reference separations are from the IERS (1996). See Table 1. 

d To facilitate comparisons between Aa and A<5, we have multiplied the Aa (ms) values by 15 • cos <5i803+784 to obtain Aa in 
mas. 
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Fig. 4. Postfit residuals of the phase delays (ns) of 1150+812 for the subset of baselines involving the following 
antennas: Effelsberg (B), Medicina (L), Fort Davis (F), and Los Alamos (X). The solution utilized corrections for 
ionosphere effects determined from our dual-frequency-band observations. For that reason, the baselines involving 
Effelsberg show a large time gap, since this antenna observed in dual-frequency mode after 23:28. Similarly, the 
Medicina antenna shows a gap until 23:48, since the phase-delay connection for this antenna appeared uncertain in 
this time range. One ambiguity interval at 8.4 GHz is 119 ps (shown at the right margin). The error bars shown are 
the standard deviations of the phase delays, scaled so that the weighted-mean-square of the postfit residuals is unity 
for each baseline (see text). 



fit residuals was nearly unity. As a result, the error bars 
assigned to the differenced phase delays were, in general, 
smaller by a factor ~ 1.2 with respect to the error bars 
assigned to the phase delays of 1150+812 (see Figs. |i]and 
||). We studied the effect of varying the relative weights 
of our two data types and verified that such variations 



did not change our result significantly. In fact, any ra- 
tio of the weights of the differenced to the undifferenccd 
data between 0.1 and 10 with respect to our "nominal" 
solution did not alter our estimate of the angular separa- 
tion significantly. Figs. ^ and |^ also illustrate a virtue of 
the difference phase-delay technique mentioned above: al- 
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though the separation of the sources is 14° 50' 21'.'147559 
and the cycle length of the observations is correspond- 
ingly large, the differenced phase delays are partially free 
from some incompletely modeled effects as attested by the 
larger residuals seen in Fig. |] than in Fig. |[ in particular 
for the baselines involving the Medicina antenna. 

The standard errors in Table || are larger than the dif- 
ferences between our estimates of Aq and A<5, and those 
provided by IERS, indicating that the a priori IERS values 
are quite accurate. This accuracy is not surprising, since 
both radio sources are dominated by emission from the 
core, and so the IERS values - mainly based on averages 
of results from observations at many epochs using group- 
delay astrometry - should be accurate, except for a small 
offset that may be contributed by persistent opacity effects 
in the core components. For sources with strong emitting 
features well separated from their cores, larger differences 
from the a priori IERS values can be expected, given the 
greater sensitivity of the "structure group delay" to com- 
ponents farther from the cores (Thomas 1980| ) . Those dif- 
ferences can be largely removed by using the more precise, 
structure-free phase delays. 



6. Sensitivity Analysis 

We carried out a sensitivity analysis to estimate the con- 
tributions of individual effects to the standard errors in 
the relative-position determination; see Table || for a list 
of these effects. We found that the primary source of error 
is due to the uncertainty in our knowledge of the neutral 
atmosphere, which by itself contributes a standard error 
of about 0.4 mas in right ascension and about 0.6 mas in 
declination to the estimate of the position of 1803+784. 
To be conservative in our estimate of the error contributed 
by the neutral atmosphere, we took for each site the larger 
of the following two values: the root-sum-square of the de- 
lay contribution from all atmosphere parameters for that 
site (one every three hours), or the algebraic sum of the 
signed partials. In this way, we made some allowance for 
the potential effects of possible correlations among the es- 
timates of the atmosphere parameters for each site. The 
second largest contribution to the total error budget of 
each coordinate comes from the corresponding standard 
error in the a priori position used for the reference source 
(1150+812). This large contribution is due to the geom- 
etry of the sources on the sky. For sources with smaller 
angular separations, the contribution of the a priori un- 
certainties of the coordinates of the reference source to 
the standard error in relative position is usually a small 
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fraction of the a priori uncertainty of the reference co- 
ordinates (e.g., Guirado et al. 1995b ), because there is 



a positive correlation between the right ascensions, and 
similarly the declinations, of the two sources. However, be- 
cause 1150+812 and 1803+784 are within 12° of the north 
celestial pole and lie almost exactly six hours apart in right 
ascension, an error in the position of one of the sources in 
right ascension (or declination) translates into a similar 
error for the other source in declination (or right ascen- 
sion). (See Table |5].) A covariance analysis confirms this 
"tradeoff" , showing a fairly strong correlation between the 
estimates of the right ascension of one source and the 
declination of the other source (p(aii5o+8i2, <5i803+784) ~ 
-0.7; jo(5ii5 0+ 8i2, "1803+784) ~ 0.7), but little correlation 
between the estimates of the two right ascensions (and 
between the two declinations). Moreover, the nearly com- 
plete cancellation of the geometric errors usually obtained 
in astrometric studies of radio sources with smaller angu- 
lar separations is not so effective in the present case. For 
example, the uncertainties in UT1-UTC values do con- 
tribute significantly to the standard error in the estimate 
of the angular separation. The remaining significant con- 
tributions to the standard error are attributable to the 
uncertainties in the estimates of the antenna locations. 

Software limitations prevented us from taking into ac- 
count in our analysis such effects as ocean and atmosphere 
loading, tidal terms in polar motion and UT1, and arbi- 
trary variations in the atmosphere zenith delays. Based 
upon a surrogate test of the sensitivity of the relative 
position of a different pair of radio sources (observed in 
an unrelated experiment) to these effects, performed us- 
ing the CALC (Ma et al. |1986[ ) and SOLVK (Herring et 
al. 1990) packages, we estimated that the known limita- 
tions of our software contribute an uncertainty of about 
0.4 mas to each coordinate of our estimate of the sources' 
relative position. 



7. Summary 

We observed the strong radio sources 1150+812 and 
1803+784 with a VLBI array on 18 November 1993, an 
epoch of mild solar activity. The antennas recorded data 
simultaneously at 8.4 and 2.3 GHz, which allowed us to 
estimate the ionosphere contributions to the phase de- 
lays; we also used TEC values from GPS measurements 
to estimate such contributions. We estimated, and thereby 
partly removed from the phase-delay data, contributions 
due to the dry and wet components of the atmosphere, 
and the brightness structure of the sources. The phase- 
connection process, required to extract the precision in- 
herent in the phase delays, did not pose special difficulties 
for our 7 min cycle time, indicating that spatial and tem- 
poral fluctuations due to the atmosphere and the iono- 
sphere were not large enough to prevent reliable phase 
connection in either frequency band. 



Table 4. Sensitivity of the estimated position of 1803+784 
to the values of other model parameters included in the 
weighted-least-squares analysis. 



Parameter 




Standard 
Deviation 3 


SAa 
( M as) b 


8AS 
(/i as) 


CW150+812 C 




174 fias 


-6 


-158 


5ll50+812 C 




140 fias 


135 


2 


Atmosphere 


r 


u.i / ns/ o nr 


— IZo 




zenith delays d : 


I 


0.17 ns/3 hr 


-150 


-87 




L 


0.17 ns/3 hr 


-267 


117 




O 


0.17 ns/3 hr 


111 


415 




X 


0.17 ns/3 hr 


-45 


227 


Site coordinates : 


B 


2 cm 


60 


40 




F 


2 cm 


30 


35 




I 


2 cm 


42 


20 




L 


2 cm 


69 


32 




O 


2 cm 


34 


76 




X 


2 cm 


24 


39 


Earth's pole c : 


X 


0.7 mas 


18 


-1 




y 


0.7 mas 


8 


21 


UT1-UTC C 




0.04 ms 


-126 


-95 


Root- Sum- Square 






412 


616 



a Each value in the SAa and 5 AS columns was obtained by 
computing the effect of altering the a priori value of each 
parameter by one standard deviation on the relative-position 
determination. 

b We have multiplied /is by 15 ■ cos <5i803+784 to obtain /ias for 
the penultimate column. 



: The standard deviations are as provided in IERS ( 1996 ). For 
the sensitivity analysis, we adopted standard deviations for 
the co ordina tes of 1150+812 twofold larger than the quoted 
IERS (|l996| ) errors. 

1 We specified the atmosphere delays at each site by a 
piecewise-linear function characterized by values at epochs 
about three hours apart. For each such epoch, we used 0.17 
ns as the standard error. The quoted sensitivity values take 
into account the total contribution, for each antenna, from all 
atmosphere parameters. Effelsberg (B) is not includ ed, since 
we solved for its atmosphere parameters. (See Sect. 4.2.) 

: As given in Table 1. The 2 cm standard deviation refers to 
each of the three coordinates for each antenna site. The 
sensitivity values account for the errors in all three 
coordinates, assumed to be independent. 



We then estimated the relative position of the sources 
(Table ||) via a weighted-least-squares analysis of a com- 
bined data set of undifferenced and differenced phase de- 
lays. The estimates resulting from use of GPS-based TEC 
values to account for ionosphere effects agree with those 
obtained from our dual-frequency-band VLBI measure- 
ments to within the standard errors for each method. 



Ros et al. (1999, 2000) also successfully used GPS-based 
TEC values to remove the ionosphere's contribution to the 
phase delays at 8.4 GHz. The checking of their removal in 
that experiment was limited to the precision of the dual- 
frequency-band (8.4 and 2.3 GHz) group delays, since Ros 
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et al. could not reliably connect their phase delays at 2.3 
GHz. In our case, we successfully connected the phase de- 
lays at both frequencies, and thus for the first time com- 
pared TEC ionosphere delays against phase-delay-based 
ionosphere corrections (Fig. ||). This agreement supports 
the use of single-band (> 8.4 GHz) VLBI observations for 
astrometry purposes, an option of particular interest for 
experiments to be carried out at epochs of strong solar ac- 
tivity, when phase connection for a frequency band < 2.3 
GHz is likely to fail. Single-frequency VLBI observations 
along with GPS-derived TEC ionosphere corrections can 
be confidently used instead. 

Our standard errors shown in Table ^ are larger than 
the difference between our estimates of relative source 
position and those provided by IERS. The high level of 
agreement is likely due to the fact that both radio sources 
are strongly core-dominated. However, were a bright emit- 
ting feature present several milliarcseconds from one of 
the cores, large offsets could exist between the true core 
position and the IERS position. With VLBI phase-delay 
astrometry, such errors can be effectively reduced using 
structure-phase corrections computed from self-calibrated 
maps derived from the same VLBI observations. More- 
over, phase-delay astrometry has the advantage of allow- 
ing us to know very accurately to which feature inside each 
source we are referring when calculating relative positions. 

In summary, we have determined with submilliarcsec- 
ond accuracy the angular separation of two radio sources 
separated by almost 15°, using phase delays from dual- 
frequency-band VLBI measurements. Since within 15° of 
any given source there are almost always two or more rea- 
sonably strong radio sources, we have thus demonstrated 
that phase delays should be usable for full-sky astrome- 
try of radio sources. We have also shown that GPS-based 
measurements can be used to obtain reliable ionosphere 
corrections to the phase delays, thus demonstrating the 
feasibility of conducting VLBI solely observations at fre- 
quencies > 8.4 GHz for astrometric purposes. Future im- 
provements in the modeling of the atmospheric delay con- 
tribution, and in the knowledge of Earth rotation and pole 
position, as well as antenna location, should result in in- 
creasingly accurate estimates of the relative positions of 
sources far apart on the sky. 

We stress that phase delays are more reliably cor- 
rected for structure phase than are group delays. There- 
fore, differenced phase-delay astrometry is better suited 
than group-delay astrometry for carrying out astromet- 
ric studies of extended radio sources. To date, however, 
group-delay astrometry has been used to establish a quasi- 
inertial celestial reference frame based on estimates of the 
positions of a number of relatively compact extragalac- 
tic radio sources from many years of regular observations. 
Our results open the avenue to an alternative, potentially 
more accurate, approach, namely that of carrying out a 
suitable series of observing sessions and using difference 



phase-delay astrometry to obtain submilliarcsecond posi- 
tions for the cores of these sources. 
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